Turbulent magnetic Prandtl numbers obtained with MHD Taylor-Couette flow 
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The stability problem of MHD Taylor-Couette flows with toroidal magnetic fields is considered in 
dependence on the magnetic Prandtl number. Only the most uniform (but not current-free) field 
with Bin = i^out has been considered. For high enough Hartmann numbers the toroidal field is 
always unstable. Rigid rotation, however, stabilizes the magnetic (kink-) instability. 

The axial current which drives the instability is reduced by the electromotive force induced by 
the instability itself. Numerical simulations are presented to probe this effect as a possibility to 
measure the turbulent conductivity in a laboratory. It is shown numerically that in a sodium 
experiment (without rotation) an eddy diffusivity 4 times the molecular diffusivity appears resulting 
in a potential difference of ^ 34 mV/m. If the cylinders are rotating then also the eddy viscosity can 
be measured. Nonlinear simulations of the instability lead to a turbulent magnetic Prandtl number 
of 2.1 for a molecular magnetic Prandtl number of 0.01. The trend goes to higher values for smaller 
Pm. 

PACS numbers: 47.20. Ft, 47.65. +a 
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INTRODUCTION 

Strong enough toroidal fields that are not current- 
free become unstable due to the Tayler instability (TI, 
Because the source of the energy is the electric 
current, these (mainly nonaxisymmetric) instabilities can 
exist even without any rotation. On the other hand, it 
is becoming increasingly clear that the stability of dif- 
ferential rotation under the presence of magnetic fields 
is one of the key problems in MHD astrophysics. There, 
however, is no laboratory experiment so far and even the 
related numerical simulations of TI are very rare |5| . 
We shall demonstrate here how the TI interacts with dif- 
ferential rotation, and how it is possible to verify the 
main results in laboratory experiments. In particular, 
the theoretical results are used to propose experiments 
for measuring the turbulent diffusivity via a Tl-induced 
reduction of the electromotive force. Such experiments 
will be of high relevance as the knowledge of the mag- 
netic turbulent diffusivity is basic for many applications 
in fluid dynamics. Very often we have only limited infor- 
mations about the magnetic diffusivity. In the laboratory 
only a very small number of experiments have been done 
(see [sl,!^)- The same is true for the eddy viscosity which 
can be measured with the same experimental device so 
that finally the turbulent magnetic Prandtl number be- 
comes known for one and the same instability. 

Consider a Taylor-Couette (TC) flow with U as the ve- 
locity, B the magnetic field, u the microscopic kinematic 
viscosity and r] the microscopic magnetic diffusivity. The 
basic state in cylindrical geometry is Ur = Uz = Br = 
Bz = and 



Let 
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bB 
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Rin and i?out are the radii of the inner and outer cylin- 
ders, Qin and flout their rotation rates and Bi^ and Bout 
are the azimuthal magnetic fields at the inner and outer 
cylinders. In particular, a field of the form hs/R is gen- 
erated by an axial current only through the inner re- 
gion R < i?in, whereas a field of the form asR is gener- 
ated by a uniform axial current through the entire region 
R < Rout including the fluid. 

The magnetic Prandtl number Pm, the Reynolds num- 
ber Re and the Hartmann number Ha, 



Pm= - 



Re=^^, 

T] V 



Ha = 



(3) 



are the basic parameters of the problem where Rq = 
\^DRiu is the unit of length with D = i^out — ^in- For 
the velocity the boundary conditions are assumed as al- 
ways no-slip {u = 0). For conducting walls the radial 
component of the field and the tangential components of 
the current must vanish so that db(f)/dR-\-b(f)/ R = bn = 
at both Riu and i^out- Here u and b are the fluctuating 
components of flow and field. 

While the linear stability code works well with small 
Pm, the minimum microscopic Pm which can be handled 
with our nonlinear code is (only) 10~^. The numerics 
cannot deal with the very small magnetic Prandtl num- 
bers of liquid metals used in the laboratory (Pm < 10~^). 
Some of our results can only be obtained by extrapolation 
methods. 



2 



THE INSTABILITY MAP 

The map for the marginal instabihty is the result of 
a linear theory for liquid sodium with Pm = 10~^. The 
linearized equations of the MHD system in a TC flow 
under the presence of a toroidal field are given elsewhere, 

i. 

Figure [T] shows the results for various values oi jiB- 
For fiB = 0.5 the toroidal field is current-free between 
the cylinders. The outer cylinder is assumed as resting 
so that for vanishing magnetic field the rotation law is 
centrifugally-unstable for Re > 68. From the given pro- 
files closest to being current-free is jub = and we do not 
find that for Ha < 200 there is any sign of destabilizing 
infiuence of the magnetic field, for neither axisymmetric 
nor nonaxisymmetric perturbations. 

For certain /j^b the m = 1 mode should be unstable 
while the m = mode should be stable The values 
/iB = 1 and jj^B = 2 in Fig. [T] are examples of this 
situation. There is always a crossover point at which 
the most unstable mode changes from m = to m = 1. 
Note also that for /i^ = 1 the critical Reynolds number 
for the m = mode steadily increases, while the m = 
1 mode is suddenly decreasing for a sufficiently strong 
magnetic field. Hence, weak fields initially can stabilize 
the fiow, and stronger fields eventually destabilize via a 
nonaxisymmetric mode. Beyond Ha =150 the fiow is 
unstable even for Re = 0. 

Except for the almost current-free profile /i^ = all 
other values share the feature that there is a critical Hart- 
mann number beyond which the basic state is unstable 
even for Re = 0. Let Ha^^^ and Ha^^^ denote these critical 
Hartmann numbers for m = and m = 1, resp. 



Flat rotation laws 

The rotation law with resting outer cylinder destabi- 
lizes the magnetic field. The question arises what hap- 
pens for those fiat rotation laws which are stable in the 
nonmagnetic regime. In Fig. [2] the marginal stability 
curves are also given for /j^q = 0.25 (Rayleigh limit), 
fiQ = 0.35, /iQ = 0.45 and fiQ = 1 (rigid rotation). One 
finds the instabilities more and more stabilized by the 
rotation. 

Note the massive quenching of the TI by rigid rotation. 
Even a rather slow rotation prevents the TI to destabilize 
the system. Rigidly rotating containers can keep much 
stronger fields as stable than without rotation. This ro- 
tational stabilization is modified for nonuniform rotation. 
At the Rayleigh limit, where Q oc R~^^ even a slow ro- 
tation destabilizes the system while it is stabilized for 
fast rotation. Generally, fast rotation stabilizes, slow ro- 
tation destabilizes. At a Hartmann number of (say) 50 
and at the Rayleigh line one finds for increasing rotation 
rate the regimes: stable, unstable, stable. The critical 



Reynolds numbers of the sequence are ~ 300 and ~ 1800 
which can easily be realized in the laboratory. A similar 
situation holds for the (quasi-Kepler) rotation law with 
/iQ = 0.35 while for rotation laws with /j^q > 0.45 only 
the rotational stabilization can be observed. 



Electric currents 

For experiments the electric currents must not be too 
strong. As an upper limit currents with 10-15 kA shall 
be considered. In order to translate the obtained criti- 
cal Hartman numbers into amplitudes of electrical cur- 
rents we apply our results to liquid sodium with a den- 
sity of 0.92 g/cm^, a microscopic magnetic diffusivity of 
810 cm^/s and a magnetic Prandtl number of 10~^. For 
gallium-indium-tin the necessary currents are stronger by 
a factor of 3.15 (see [10]). 

Let /axis be the axial current inside the inner cylinder 
and /fluid the axial current through the ffuid (i.e. be- 
tween inner and outer cylinder). Then the toroidal field 
amplitudes at the inner and outer cylinders are 

-^axis 7-, (-^axis H~ -^fluid) 



Si, 



Bo 



(4) 



5i?out 

measured in cm. Gauss and Ampere. Expressing /axis 
and /fluid in terms of our dimensionless parameters one 
finds 



/a: 



5Ha 



1 — 77 ' 



J _ f^B -r] J 

-'fluid — ^ -'axis- 



Table [J gives the electric currents needed to reach the less 
of Ha^^^ and Ha^^^ for fj = 0.5, and /j^b ranging from —2 
to 2 in each case. Note that for large l/i^l the current 
^fluid approaches a constant value. 



TABLE L Characteristic Hartmann numbers and electric cur- 
rents for a sodium-container (77 = 0.5) with conducting walls. 
The experiment with the almost uniform field fiB — 1 is indi- 
cated in bold. 



fiB Ha(°) Ha(^) /axis [kA] h 


uid [kA] 


-2 19.8 24.8 


0.807 


-4.04 


-1 59.3 63.7 


2.42 


-7.25 


1 00 151 


6.16 


6.16 


2 00 35.3 


1.44 


4.32 



The most interesting experiment is that with the al- 
most uniform field /i^ = 1 • For a container with a gap of 
fj = 0.5, parallel currents of 6.16 kA are necessary along 
the axis and through the fiuid. The experiment does not 
possess the weakest electric currents but both the cur- 
rents are parallel and have the same amplitudes. Figure 
[1] (middle) shows that in this case a crossing point M 
exists where the axisymmetric mode has the same char- 
acteristic Reynolds number and Hartmann number as the 
nonaxisymmetric mode with m = 1. 



3 




FIG. 1: The marginal stability curves for m = (dashed) and m ■ 
Conducting walls. The hatched domain is stable. 



1 (solid). Pm= 10 ^ ^ fj — 0.5, im — 0, and /is as indicated. 
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FIG. 2: fjiB — ^- The same magnetic constellation as in Fig. 
[T] (middle) but for fin as indicated, /in = (resting outer 
cylinder), fin = 0.25 (Rayleigh limit), /in = 0.35, /in = 0.45, 
/in = 1 (rigid rotation) . The dotted line gives the eigenvalues 
for m = which for those parameters only exists for /in = 0. 



for the eddy diffusivity, [ll'] . The magnetic fluctuations 
do not contribute to the magnetic diffusivity. This basic 
difference between both the diffusion coefficients is not 
yet proven by an experiment. The results (|6]) and ([7j) 
suggest that in turbulent magnetic fluids the effective 
magnetic Prandtl number exceeds the value 0.4 which 
was conflrmed by numerical simulations for driven MHD 
turbulence with Pm of order unity, [12]. The knowledge 
of the turbulent magnetic Prandtl number is of extraor- 
dinary meaning in fluid mechanics and geo/ astrophysics. 
For its calculation one has to measure both quantities 
simultaneously in one and the same experiment. 

Simplifying, the nonaxisymmetric components of flow 
and fleld may be used in the following as the 'fluctuations' 
while the axisymmetric components are considered as the 
mean quantities. Then the averaging procedure is simply 
the integration over the azimuth (j). It is standard to ex- 
press the turbulence-induced electromotive force (EMF) 



as 



THE EDDY DIFFUSIVITY 

We now turn to the mean-fleld concept turbulent fluids 
of electrically conducting material. It is known that the 
existence of turbulence in the fluid reduces the electric 
conductivity or - with other words - the fluctuations en- 
hance the magnetic diffusivity called the turbulent mag- 
netic diffusivity. In MHD the turbulent diffusivity is a 
much more simple quantity than the corresponding eddy 
viscosity. While the latter is also formed by the existing 
magnetic fluctuations this is not the case for the turbu- 
lent diffusivity. In a simplifled ('SOCA') approximation 
for a turbulence fleld with a correlation time Tcorr results 



S = {u X b) = —tjtcutIB 



(8) 
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for the eddy viscosity but only 



(6) 



(7) 



with the (scalar) eddy diffusivity r]T which must be pos- 
itive. In cylindric geometry the mean current cur IS has 
only a 2;-component. Hence, 



Take from Table [J that for /i^ = 1 the current through 
the fluid is positive then for negative Ez the tjt results as 
positive. This is indeed the case. We have shown that TI 
indeed provides reasonable expressions for the turbulent 
diffusivity in rotating containers, [13|]. 



Nonlinear simulations 

The absolute values for £z can only be computed with 
nonlinear simulations. The minimum possible magnetic 
Prandtl number for the code yielding robust results is 
10~^. Here the results without and with rotation but 
only for Ha = 200 are reported. The used MHD Fourier 
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spectral element code has been described earlier in more 
detail, [11, [H. 

Either M = 8orM=16 Fourier modes are used, 
two or three elements in radius and twelve or eighteen 
elements in axial direction, resp. The polynomial order 
is varied between N = S and TV = 16. 

Figure [3] shows the negative Tl-induced EMF for mag- 
netic Prandtl numbers varied between 0.01 to 1. The 
main results are that i) the EMF is always negative (t^t 
positive!) and ii) it runs with £^/Pm with the factor 
^ :^ 3 taken from the plot. The resulting EMF in physi- 
cal units is —tjEBq/D. Hence, 



V 



1.5E 4.5. 



(9) 



For the voltage difference 5U due to this EMF one finds 
6U = tjEBqH/D with H as the container height and 



6U 



D 



(10) 



with r = H/D the aspect ratio of the container. For 
= 10 cm and H = 100 cm we find for sodium 
{-^/JlopUfj :^ 8.15) the maximum value of 34 mV as the po- 
tential difference from endplate to endplate. This value 
can only be considered as an estimate basing on the scal- 
ing with 1/Pm suggested by Fig. [3l But even in the case 
that the slope of the curve decreases for smaller Pm the 
effect should be observable in the laboratory. A Hart- 
mann number of 200 requires 163 G at the inner cylinder 
(^in = 10 cm) which can be produced with an axial cur- 
rent of 8.15 kA for i^out = 2i?in. 
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FIG. 3: Maximum values of the normalized axial EMF vs 
Pm for Ha = 200 and /is = 1- The lower curve is without 
rotation, the upper curve is with rotation and shear (Re = 
500, /XQ = 0.35). 

In order to study the rotational influence also the Pm- 
dependence of the EMF under the presence of a differen- 
tial rotation is given in Fig. [3l It is ji^t = 0.35 (quasikep- 
lerian) and the Reynolds number is Re = 500. Note that 
the influence of the rotation for the Tl-induced EMF is 
surprisingly weak; with rotation the values are slightly 
higher than without rotation. 



The magnetic Reynolds number, Rm = MAX(2x)Z)/?7, 
of the fluctuations is considered next. With Fig. [H a 
rather weak magnetic Prandtl number dependence of Rm 
is found. Extrapolating the results to Pm = 10~^ gives 
in both cases a value of Rm 2.6. The associated ve- 
locity fluctuations for sodium are about 15 m/s in a gap 
of 1 cm and 1.5 m/s in a gap of 10 cm. The values are 
rather similar to those of the Riga 'a-yashchik' experi- 
ment, [16]. Even with resting cylinders it is possible to 
produce rather high (azimuthal) velocities in TI experi- 
ments. 
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FIG. 4: The same as in Fig. [3] but for the magnetic Reynolds 
number of the fluctuations. 



THE EDDY VISCOSITY 

Experiments with Tayler instability under the presence 
of differential rotation can also provide eddy viscosity 
measurements due to the angular momentum transport 
by both Reynolds stress and Maxwell stress. Within the 
diffusion approximation it is 



(11) 



for the torque in the fluid. The fluctuations of flow and 
fleld can be calculated with the code. The patterns for 
the instability-induced diffusivity values for Pm = 0.01 
are shown in Fig. [5] for the same model as used in Fig 
m i.e. i^n = 0.35, Ha = 200 and Re = 500. One flnds 
the turbulence-originated increase of the eddy viscosity 
ut/^ much larger than for the turbulent diffusivity. The 
turbulent magnetic Prandtl number Pmx = i^t/vt be- 
comes about 2.05. Similar calculations for Pm = 0.1 
lead to Pmx = 0.71 while for Pm = 1 the smaller value 
0.65 results. The results only weakly depend on the av- 
eraging procedure. The given numbers follow after aver- 
aging over the whole cylinder. The turbulent magnetic 
Prandtl number slightly increases with decreasing micro- 
scopic Pm; and for small Pm it reaches values larger than 
unity. Note the differing results of simulations with Pm 
much smaller than unity and those with Pm < 1, [12]. 
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FIG. 5: Ha = 200, Re = 500. Simulations for a flow with Pm 
resulting magnetic Prandtl number vt/tit is about 2.1. 



Small Pm are shown to produce large turbulent values 
Pmx. We cannot provide results for Pm smaller than 
0.01 so far. Only laboratory experiments utilizing Tayler 
instability in liquid metals with their very small magnetic 
Prandtl number are able to show whether the obtained 
trend is a general one. 
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